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Abstract 

We introduce a transfer matrix formalism for the (annealed) Ising 
model coupled to two-dimensional causal dynamical triangulations. 
Using the Krein-Rutman theory of positivity preserving operators we 
study several properties of the emerging transfer matrix. In partic- 
ular, we determine regions in the quadrant of parameters /3, > 
where the infinite- volume free energy converges yielding results on the 
convergence and asymptotic properties of the partition function and 
the Gibbs measure. 

2000 MSG. 60F05, 60J60, 60J80. 

Keywords: causal dynamical triangulation (CDT), Ising model, par- 
tition function, Gibbs measure, transfer matrix, Krein-Rutman theory 



1 



1 Introduction. A review of related results 



In the study of two-dimensional quantum gravity and non-critical string the- 
ory, models of discrete random surfaces play an essential role. 

In the 1980s, so-called dynamical triangulations (DT) were introduced to 
define a Euclidean path integral for two-dimensional quantum gravity (see [1] 
for an overview). In particular, the partition function has been determined 
as a sum over all possible triangulations of a sphere where each configuration 
is weighted by a Boltzmann factor e"^'^', with |T| standing for the size of 
the triangulation and ji being the cosmological constant. The evaluation of 
the partition function was reduced to a purely combinatorial problem that 
can be solved with the help of the early work of Tutte [21 [3] ; alternatively, 
more powerful techniques were proposed, based on random matrix models 
(see, e.g., [4J) and bijections to well-labelled trees (see E])- One can 
then pass to a continuum limit by taking the number of triangles to infinity. 
An interesting property of the resulting "quantum geometry" is its fractal 
structure as illustrated in Figure [i] (a). In the physical literature such fractal 
structures are called "baby universes", and they completely dominate the 
continuum limit leading to a fractal dimension = 4. 

From a probabilistic point of view there has recently been an increasing 
interest in DT, most notably through the work of Angel and Schramm on a 
uniform measure on infinite planar triangulations [7j, as well as through the 
work of Le Gall, Miermont and collaborators on Brownian maps (see [8j for 
a recent review) . 

From a physical point of view it is interesting to study various models 
of matter, such as the Ising model, coupled to the DT. The calculation of 
the partition function in this case also reduces to a combinatorial problem. 
It was first solved in O [lO] by using random matrix models and later by 
using a bijection to well-labelled trees [H]. It is interesting that the solution 
here is much simpler than in the case of a fiat triangular or square lattice as 
given by Onsager [12j. Further, one can see that the critical exponents in the 
case where the model is coupled to DT differ from the Onsager values. This 
is related to the strong back-reaction of the Ising model with the quantum 
geometry. In particular, the spin clusters energetically prefer to sit within 
baby universes since those are connected to the main universe through a very 
short so-called bottleneck boundary (see Figure [i] (b)). The spins increase 
the fractal structure leading to a change in the values of the critical exponents 
at the critical temperature. 
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(a) 



(b) 



Figure 1: (a) A section of a typical planar triangulation of the DT ensemble 
illustrating the fractal structure of the quantum geometry, (b) Illustration 
of the baby universes and the formation of spin clusters within them. Each 
baby universe corresponds to a fractal structure as in part (a) drawn out of 
the plane. 

In a continuum framework one attempts to understand the resulting the- 
ory as Liouville theory coupled to a conformal field theory with central charge 
c = 1/2. Furthermore, this leads to a simple algebraic identity (the KPZ- 
relation) between the critical exponents of the Ising model on a fiat lattice 
and the critical exponents of the Ising model coupled to DT [13j. 

While DT has a very rich mathematical structure which very recently has 
been related to the SLE (the Schramm-Loewner evolution) and level curves 
of a Gaussian free field [14j, from the point of view of quantum gravity its 
fractal structure leads to causality-violating geometries that are arguably 
non-physical. This led to the development of so-called causal dynamical 
triangulations (CDT) by Ambj0rn and Loll [15j, to define the Lorentzian 
gravitational path integral. A causal triangulation is formed by triangula- 
tions of spatial strips as illustrated in Figure [2} Note that the left and right 
boundaries of the spatial strip are periodically identified. 

A first analytical solution of two-dimensional (pure) CDT was obtained 
in [15] where it was shown that the resulting quantum geometry, while still 
random, is much more regular than in the case of a DT, leading to a fractal 
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Figure 2: A strip of a causal triangulation of 5 x [j^j + 1]. 



dimension of d = 2. From a probabilistic point of view, we would like to 
note that a uniform measure on infinite causal triangulations UICT has been 
recently introduced in [IS^^TT]. 

An interesting question is: What are the properties of the Ising model 
coupled to a CDT ensemble? As was said above, the CDT ensemble is more 
regular than that of the DT, but it is still random and allows for a back- 
reaction of the spin system with the quantum geometry. Monte Carlo simu- 
lations [18j (see also [tl9||20j) give a strong evidence that critical exponents 
of the Ising model coupled to CDT are identical to the Onsager values. 

While recently much progress has been made in the development of an- 
alytical techniques for CDT particularly random matrix models 
[231 ES EH], and their application to multi-critical CDT [26l Ell EH], the 
causality constraints still makes it difficult to find an analytical solution of 
the Ising model coupled to CDT. For the quenched Ising model coupled to 
two-dimensional CDT some progress has been made in proving the existence 
of a phase transition [29j. 

In this article we develop a transfer matrix formalism for the annealed 
Ising model coupled to two-dimensional CDT. Spectral properties of the 
transfer matrix are rigorously analysed by using the Krein-Rutman theo- 
rem [30] on operators preserving the cone of positive functions. This yields 
results on convergence and asymptotic properties of the partition function 
and the Gibbs measure and allows us to determine regions in the parameter 
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quarter-plane where the partition function converges. 

Outline: The article is organised as follows. Section [2] contains basic 
definitions. Section |3] gives a summary of the transfer-matrix formalism for 
CDTs. In Section |4] we introduce the transfer matrix for the Ising model 
coupled with a CDT and state and comment on our main results. This is 
followed with concluding remarks in Section [5} An Appendix contains the 
proofs of the stated results. 

2 Definitions 

We will work with rooted causal dynamic triangulations of the cylinder Cn = 
5 X [0, A"], = 1, 2, ... , which have A bonds (strips) S x + 1]. Here 
S stands for a unit circle. The definition of a causal triangulation starts by 
considering a connected graph G embedded in Cn with the property that all 
faces of G are triangles (using the convention that an edge incident to the 
same face on two sides counts twice, see [17j for more details). A triangulation 
t of Cat is a pair formed by a graph G with the above propetry and the set 
F of all its (triangular) faces: t = (G, F). 

Definition 2.1. A triangulation t ofC^ is called a causal triangulation (CT) 
if the following conditions hold: 

• each triangular face of t belongs to some strip S x [j^j + 1]; j = 
1, . . . , A — 1^ and has all vertices and exactly one edge on the boundary 
{S X {j}) U{Sx {j + 1}) of the strip S x + 1]; 

• if kj = kj{t) is the number of edges on S x {j} , then we have < kj < 
oc for a// j = 0, 1, . . . , A — 1. 

Definition 2.2. A triangulation t of Cn is called rooted if it has a root. The 

root in the triangulation t is represented by a triangular face t of called 
the root triangle, with an anticlock-wise ordering on its vertices (x, z) where 
X and y belong to x {0}. The vertex x is identified as the root vertex and 
the (directed) edge from x to y as the root edge. 

Definition 2.3. Two causal rooted triangulations ofCN, say t = (G, F) and 

X! = {G\ F'), are equivalent if there exists a self-homeomorphism ofCN which 
(i) transforms each slice aS-^ x {j}^ j = 0, . . . , A — 1 to itself and preserves its 
direction, (ii) induces an isomorphism of the graphs G and G' and a bisection 
between F and F' , and (Hi) takes the root of t to the root oft/ . 



5 




where t{i) is a causal triangulation of the strip S x[i^i + l]. The latter means 
that each t{i) is described by a partition of 5 x [i, i + 1] into triangles where 
each triangle has one vertex on one of the slices S x {i}, S x {i + 1} and two 
on the other, together with the edge joining these two vertices. The property 
of consistency means that each pair (t(i),t(i + 1)) is consistent, i.e., every 
side of a triangle from t{i) lying in S x {i + 1} serves as a side of a triangle 
from t(i + 1), and vice versa. 

The triangles forming the causal triangulation t(i) are denoted by j), 
1 < J < '^(t(O) where, n{t{i)) stands for the number of triangles in the 
triangulation t{i). The enumeration of these triangles starts with what we 
call the root triangle in t(i); it is determined recursively as follows. First, we 
have the root triangle t(0, 1) in t(0) (see Definition 2.2). Take the vertex of 
the triangle t(0, 1) which lies on the slice 5 x {1} and denote it by x^ This 
vertex is declared the root vertex for t(l). Next, the root edge for t(l) is the 
one incident to and lying on 5 x {1}, so that if is its other end and is the 
third vertex of the corresponding triangle then x\y\ z' lists the three vertices 
anticlock-wise. Accordingly, the triangle with the vertices x'^y\z' is called 
the root triangle for t(l). This construction can be iterated, determining the 
root vertices, root edges and root triangles for t(i), < i < — 1. 

It is convenient to introduce the notion of "up" and "down" triangles (see 
Figure [2]). We call a triangle t G t(i) an up-triangle if it has an edge on the 
slice S X {%] and a down-triangle if it has an edge on the slice S x By 
Definition 2.1, every triangle is either of type up or down. Let n^^{t{i)) and 
^do(t(i)) stand for the number of up- and down-triangles in the triangulation 



Note that for any edge lying on the slice S x {i] belongs to exactly two 
triangles: one up-triangle from t{i) and one down-triangle from t(i — 1). This 
provides the following relation: the number of triangles in the triangulation 
t is twice the total number of edges on the slices. More precisely, let be 
the number of edges on slice S x {%]. Then, for any i = 0, 1,...,A— 1, 



n{t{i)) = n„p(t(i)) + ridoitii)) = rC + n' 



(2.1) 



implying that 



N-1 



N-1 




(2.2) 
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There is another useful property regarding the counting of triangulations. 
Let us fix the number of edges and n^+^ in the shces S x {i} and Sx{i + 1}. 
The number of possible rooted CTs of the slice S x + with up- and 
n^+^ down-triangles is equal to 



n^ + n^+i - 1\ _ / n(t(i)) - 1 



(2.3) 



3 Transfer matrix formalism for pure CDTs 

We begin by discussing the case of pure causal dynamical triangulations, as 
was first introduced in (see also [31 J for a mathematically more rigorous 
account). 

The partition function for rooted CTs in the cylinder Cn with periodical 
spatial boundary conditions (where t(0) is consistent with t(A^ — 1)) and for 
the value of the cosmological constant /x is given by 

N-l 

Z^{^) = 5]e-^'^® = exp{-^5]n(t(.))}. (3.1) 

t (t(0),...,t(Ar-l)) i=0 



Using the properties (2.2) and (2.3) we can represent the partition function 



(3.1) in the following way 

N-l N-l 



z.M= E -p{-2fE"in("'trrr')- (32) 

Moreover, Zn[ii) admits a trace-related representation 

Z^(/x)=tr ([/^). (3.3) 

This gives rise to a transfer matrix U = {u{n^n^)}n,n'=i,2,... describing the 
transition from one spatial strip to the next one. It is an infinite matrix with 
strictly positive entries 

«(n,n')=('' + ;;^^^)^"+"'. (3.4) 

For notational convenience we use the parameter g = (a single-triangle 
fugacity). The entry u{n^n') yields the number of possible triangluations 
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of a single strip (say, S x [0, 1]) with n lower boundary edges (on S x {0}) 
and upper boundary edges (on S x {1}). See Figure [2| The asymmetry 
in n and is due to the fact that the lower boundary is marked while the 
upper one is not. However, a symmetric transfer matrix U = {u{n^n')} 
can be introduced, associated with a strip where both boundaries are kept 
unmarked: 

S(n, n^) = n~^u{n^ n'). (3.5) 

The A^-strip Gibbs distribution Vn assigns the following probabilities 
to strings (n^, . . . ^n^~^) with the number of triangles > 1 for all i = 
0,...,iV-l: 

N-1 N-1 / i, 

^"<" z^-" -2"^"* n(" ;r_r ■ <^-«) 

^^^^^ i=0 i=0 ^ / 

We state two lemmas featuring properties of matrix U : 

Lemma 3.1. For any g > the matrix U and its transpose have an 
eigenvalue A = A{g) given by 

A(^?)= [(l-yr^)/(2<?)]'. (3.7) 

The corresponding eigenvectors 

(f) = {(t){n)}n=i,2,... and 0* = {(/>*(n)}^^i,2,... 

have entries 

<P{n) = n{A{g)y, {n) = {K{g)r • (3.8) 
Proof. A direct verification shows that 

^u{n,n')n'K''\g)=nK''+\g) and ^K''{g)u{n,n') = K'''+\g). 

n' n 

(In fact, each of these relations implies the other.) See Theorem 1 in [3T]. □ 
Lemma 3.2. For any fixed n and any g < I (equivalently, fi > 0) one has 

Y,u{n,n') = {-^J{1 - (1 - 9T)- (3.9) 
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Proof. The proof again follows from a straightforward verification. □ 

A transfer-matrix formalism of Statistical Mechanics predicts that, as 
^ oc, the partition function is governed by the largest eigenvalue A of 
the transfer matrix: 

Z^(^) = tr [/^ - A^ (3.10) 
We make this statement more precise in the statements of Lemma |3.3| and 



Theorem 3T below. Here the symbol i'^ stands for the Hilbert space of square- 
summable complex sequences (infinite-dimensional vectors) ip = {il^{n)}n=i,2,... 
equipped with the standard scalar product {ip'^ip") = Yln^'i^)^ (^)- 
cordingly, the matrices U and are treated as operators in l^. 

Lemma 3.3. For any g < 1/2 (equivalently fi > ln2) the following state- 
ments hold true: 

1. U and are bounded operators in P preserving the cone of positive 
vectors; 

2. The sum ^/ '^(^, '^0 ^ Consequently, U and have 

tr {UU^) = tr {U^U) < oc, 

i.e., U and are Hilbert- Schmidt operators. Therefore, W N > 2, 
and {U^^ are trace-class operators. 

3. The maximal eigenvalue A = K{g) of U in P is positive, coincides 
with the maximal eigenvalue of and is given by Eqn (3.1). The 
corresponding eigenvectors 0, 0* G are unique up to multiplication 
by a constant factor and given in Eqn (3.8). 

4. The following asymptotical formulas hold as N ^ 00: 

^tr(n. ^tr((C/-n^l, 

and, V vectors ip' ^ip" G , 

^(V^',t/V) = (V'',0)(r,V'"), 

where the eigenvectors (f> and (f>* are normalized so that {4>,4>*) = 1. 
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Theorem 3.1. For any g < 1/2 the following relation holds true: 



lim —log ZN{g) = log A 



(3.11) 



with A = A{g) given in (3.1). Further, the N -strip Gibbs measure Vn 
converges weakly to a limiting measure V which is represented by a posi- 
tive recurrent Markov chain on Z+ = {1, 2, . . with the transition matrix 
P = {P{n^n')}n=i^2,... cind the invariant distribution n. Here 



P{n,n') = 



and 



u{n^ n')(l){n') 
A(l){n) 

(j)''{n)(f){n) 



where ((){n) and (l)*{n) are as in ^3^. 



Proof. The proof is a consequence of Lemma |3.1| and |3.3| and the Krein- 
Rutman theory [30]. □ 



4 An Ising model coupled to the CDT: state- 
ment of results 

4.1 The model 

With any triangle from a triangulation t we associate a spin taking values 
±1. An A^-strip configuration of spins is represented by a collection 

g:={a{0),a{l),...,(T{N-l)) 

where a{i) = a{t{i)) is a configuration of spins cr{i^j) over triangles t{i^j) 
forming a triangulation t(i), 1 < j < n{t{i)). We will say that a single- 
strip configuration of spins cr{i) is supported by a triangulation t{i) of strip 

5 X [i^i + 1]. We consider a usual (ferromagnetic) Ising-type energy where 
two spins cr{i^j) and cr{i'^j') interact if their supporting triangles t{i^j)^ 
t{i' ^f) share a common edge; such triangles are called nearest neighbors, and 
this property is refiected in the notation ((7(i, j), cr(i', /)), where we require 
{)<i<i'<N — 1. Thus, in our model each spin has three neighbors. 
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Moreover, a pair {a{i^ j)^a{i' ^ f)) can only occur for — z < 1 or z = 0, 
i' — N — 1. Formally, the Hamiltonian of the model reads: 

H(o:) = - <^.3W,f)- (4.1) 

{a{i,j),a{i>,j')) 

We will use the following decomposition: 

H(o:) = ^(^(^)) + E ^(^W' + 1))' (4-2) 

where we assume that cr(0) = cr{N) (the periodic spatial boundary condi- 
tion). Here H{a{i)) represents the energy of the configuration cr(i): 

H{a{i))^- Y ^(^'iM^'/)- (4-3) 

{a(i,j),a(i,j')) 

Further, V{a{i),a{i + 1)) is the energy of interaction between neighboring 
triangles belonging to the adjacent strips S x [z, i + 1] and S x [z + 1, z + 2]: 

y(or(z), or(z + 1)) = - ^('^JM' + 1^/)- (4.4) 

The partition function for the (annealed) A^-strip Ising model coupled to 
CDT, at the inverse temperature ^ > and for the cosmological constant /i, 
is given by 

N-l 

S;v(/^,/?)= J2 exp{-/. J]n(t(z))} (4.5) 

(t(0),...,t(Ar-l)) i=0 

N-l 

X J2 n - PViaii), <T{i + 1))}. 

{cr{o),...,cr{N-i)) i=o 

Here n(t(z)) stands for the number of triangles in the triangulation t(z). Like 
before, the formula 

Hjv(/^,/5)=trK^ (4.6) 

gives rise to a transfer matrix K with entries i^((t,cr), (t',cr')) labelled by 
pairs (t, cr), (t', cr') representing triangulations of a single strip (say, Sx [0, 1]) 
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and their supported spin configurations which are positioned next to each 
other. FormaUy, 

ir((t,o-),(t',cr')) = lt.t'exp{-|(n(t) + n(t'))} (4.7) 
X e^p{-^{HicT) + Hia'))-/3Via,a')]. 

As earlier, n(t) and n(t') are the numbers of triangles in the triangulations 
t and t^ The indicator It^t^ means that the triangulations t,t' have to be 
consistent with each other in the above sense: the number of down-triangles 
in t should equal the number of up-triangles in t', and an upper-marked edge 
in t should coincide with a lower-marked edge in triangulation t^ It means 
that the pair (t, t') forms a CDT for the strip S x [0, 2]. 

We would like to stress that the trace tr in ( |4.6[ ) is understood as 
the matrix trace^ i.e., as the sum Yl,to'-^^^\i^^^)^i^'>^)) diagonal 
entries K(^^((t,cr), (t,cr)) of the matrix K^. (Indeed, in what follows, the 
notation "tr" is used for the matrix trace only.) Our aim will be to verify 



that the matrix trace in (4.6) can be replaced with an operator trace invoking 
the eigenvalues of K in a suitable linear space. 

As before, we can introduce the A^-strip Gibbs probability distribution 



associated with formula (4.5): 



P^^((t(0), o-(O)), . . . , (t(iV - 1), ct{N - 1))) (4.8) 

N-l 

= n exp{-^n(t(^)) - /3/^(a(^)) - PV{cr{i), a{i + 1))}. 

Consider several special cases of interest. 

The case /3 0. This is the first term of the so-called high temperature 
expansion fT8]. Here one has 



S(^,0) = exp{-/i5]n(t(^))} 1 

(t(o),...,t(Ar-i)) i=o (cr(o),...,cr(Ar-i)) 

N-l N-l / 

Y exp{-2(^-ln2)J]n^}n _ , 



+ rf^^ — 1 



= Z^(/i-ln2); cf. (3.1). 
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The condition /i — ln2 > In 2 which guarantees properties hsted in 
Lemma 13.31 and Theorem 13.11 resuls in 



/i > 2 In 2. 



(4.9) 



Thus, Eqn. (4.9) yields a sub-criticality condition when /3 = 0. 



The case /3 oc. Observe that for any triangulation t = (t(0), . . . , t(A^ — 
1)) there are two ground states: aU spins +1 and aU spins —1, with 
the overaU energy equals minus three half times the total number of 
triangles: —3/2 X^^q^ n(t(i)). Discarding all other spin configurations, 
we obtain that 

S(/i,^) > 

where 



J2 2exp{{-^,+ -^)^n{t(^))} 

t(0),...,t(A^-l) i=0 

2 Y: exp{-2(^-|/3)j:n^} 



nO>l,...,n^-l>l 

3 



n' + n*"*"-^ — 1 
n' -1 



= 2Z^(^--/5 



where exp 



is the energy of the (+)-configuration (or. 



equivalently, the (— )-configuration). For /3 large, we can expect that 
/3) ^ /3). Then the critical inequality 



yields 



/x--/5>ln2 



/X > ln2 + 



(4.10) 



Equation (4.10) gives a necessary (and probably tight) criticality con- 



dition for the Ising model under consideration for large values of A 
similar result was obtained in [18j . 
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The case < ;5 < oo. Firstly, we note that for any fixed triangulation t the 
energy of any spin configuration a on t will be bigger or equal than 
the energy of a pure configuration (all +s or all — s): 

j 3 

3 ^"^ 
> — -#(of all triangles in t) = — 3 

where is the number of edges in the ith level S x {i}, i = 0, 1 . . . , — 
1. Thus, for any /5 > the inequality < holds true, 

where 

3 



S*(/i,/5) = exp{(-^ + -/5 + ln2) J]n(t(^))} 

(t(0),...,t(iV-l) i=0 

exp{-2(^--/3-ln2)5]n*} 



= Z^,(^-^/3-ln2) . 

Hence, the inequality 

/i--/3-ln2>ln2 or jU>21n2 + -/5 (4.11) 

provides a sufficient condition for subcriticality of the Ising model under 
consideration. 

4.2 The transfer-matrix K and its powers 

The main results of this article are summarized in Lemma |4T] and Theorems 



4.1 and 4.2 below. Let us start with a statement (see Proposition 4.1 below) 
which merely re-phrases standard definitions and explains our interest in the 
matrices K, K^, K^K, KK^ and their powers. Cf. Definition 2.2.2 on p. 83, 
Definition 2.4.1 on p. 101, Lemma 2.3.1 on p. 85 and Theorem 3.3.13 on p. 139 
in [32J). 
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We treat the transfer-matrix K and its transpose as linear operators 
in the Hilbert space i^_c (the subscript T-C refers to triangulations and spin- 
configurations) . The space ^x-c formed by functions i/j = {'?/^(t,cr)} with 
the argument (t,cr) running over single-strip triangulations and supported 
configurations of spins, with the scalar product '0'')t_c = St a ^)^"{^-> ^) 
and the induced norm ||'0||t-c- The action of K in ^x-c^ basis formed 

by Dirac's delta- vectors 5(t,cr)5 is determined by 

(KV')(t,o-)= J]i^((t,o-),(t',o-0)V'(t',o-'); (4.12) 
t'.cr' 

in following we use the notation K, K^, etc., for the matrices and the cor- 
responding operators in ^t-c- Accordingly, the symbols ||K||t-C5 ||K^||t-c 
etc. refer to norms in ^\_q. 

Given n = 1,2, . . ., suppose that the operator (respectively, (K^)^) 
is of trace class . Then the following series absolutely converges: 



Af ^ ( respectively, A^f^ ) , (4.13) 



where A^-^^ (A*^-^^) runs through the eigenvalues of ((K^)^), counted 



with their multiplicities. In this case the sum (4.13) is called the operator 
trace of (respectively, (K^)^) in ^\_q. We adopt an agreement that 
the eigenvalues in ( |4.13D are listed in the decreasing order of their moduli, 
beginning with Aq^^ (A*o^^) . 

Set |K-| = V(KT)"K- and \ {K^Y\ = ^K- (K^)". 

Proposition 4.1. For any positve integer r, the following inequalities are 
equivalent: 

tr((K^(KT)T=tr((KT)^KO<oo and 

tr|K2-| =tr|(KT)2-| < oc. ^ ' 



Moreover, each of the inequalities in (4.14) implies thatM N > 2r, the opera- 
tors and (K^)^ are of trace class in^\_Q. Hence, for N > 2r, the matrix 
traces tr (K^) and tr((K^)^) are finite and coincide with the corresponding 
operator traces in ^x-c- 



Theorem 4.1. Suppose that the condition (4-14) satisfied with r = 1. 



Then the following properties of transfer matrix K are fullfilled. 
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1. The square K and its transpose (K ) are trace-class operators in 

2. K and have a common eigenvalue, A = Ao(/3,/x) > such that 
the norms ||K||t-c = ||K^||t-c = A. Furthermore, and (K^)^ 
have the common eigenvalue A = Aq = A*o such that the norms 
||K2||t_c = ||(KT)2||t_c = A2 . 

3. K is a simple eigenvalue ofK. and K.^ , i.e., the corresponding eigen- 
vectors = {0(t,cr)} and 0* = {0*(t,cr)} are unique up to multi- 
plicative constants. Moreover, and 0^ can he made strictly positive: 
0(t, o"), 0^(t, cr) > V (t,cr). Furthermore, A is separated from the 
remaining singular values and the remaining eigenvalues of K and 

hy a positive gap. The same is true for A^ and and (K^) . 



Proof of Theorem 4.1. Because the entries X((t, cr), (t^, cr')) are non- 
negative, the condition ( |4.14D with r = 1 means that 



^ X2((t,o^),(t^o^0)<oo, (4.15) 

(t,o■),(t^o■o 

that is, K and are Hilbert-Schmidt operators. It means that the operator 
KK^ has an orthonormal basis of eigenvectors and the series of squares of 
its eigenvalues (counted with multiphcities) converges and gives the trace 
trT-c(KK^). Consequently, the operators K and are bounded (and even 
completely bounded) and and (K^)^ are of trace class. The latter fact 
means that the matrix trace of the operator coincides with its operator 
trace in and the same is true of (K^)^. In addition, the operator has 

the property that its matrix entries /^(^^((t, cr), (t^, cr')) are strictly positive: 

x(2)((t,a),(t^aO) = 5] i^^((t,a),(t,5^))x((t,5^),(t^aO (4.16) 

(t,5^) 

The Krein-Rutman theory (see ^0], Proposition YIV) guarantees that both 
K and have a maximal eigenvalue A that is positive and non-degenerate, 
or simple. That is, the eigenvector of K and the eigenvector 0* of 
corresponding with A are unique up to multiplication by a constant, and all 
entries 0(t,cr) and 0*(t,cr) are non-zero and have the same sign. In other 
words, the entries 0(t,cr) and 0*(t,cr) can be made positive. The spectral 
gaps are also consequences of the above properties. □ 
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Set: 

A(/i, /5) = (m^ + 1) (cosh 2/5) j 1 
where c and m are determined by 




1)2 



(cosh 2/3)2 (m2 + l)2 



(4.17) 



c = 



exp(/3 — yu) 



e2/3(l-exp(^-^))2-e-2A' 
e2^ + (l-e^^)exp(-(/5 + ^)). 



(4.18) 
(4.19) 

(4.20) 

(4.21) 

implying the assertions of Proposition j| and Theorem j[ Moreover^ the 
condition (4.14) implies (4.20) 



Lemma 4.1. For any > 5t^c/i i/^ai 

A(/x,/3) < 1, 



the condition (4.14) Z5 satisfied for r = 1: 

tr(KK^) = tr(K^K) < oo and tT\K^\ = tT\{K^f\ < oc, 



The proof of Lemma |4.1| is given in the next section. Here we only remark 
that the proof is based on the foUowing representation of the trace (4.21): 
there exists matrices Q^Qm (see formulas ( A. 11]) and ( A.12D ) with positive 
entries and of size 4x4 such that (see formula (A. 15)) 



tr(KK^ 



tr 



((E«')ft 



+ 



k>l 



The convergence of the matrix series X]iEc>i is equivalent to the condition 
that the maximal eigenvalue of the matrix Q is less then 1. This is exactly 
the condition (4.20). 



Theorem 4.2. Under condition (4-20), the following limit holds: 



lim ^logS^(/3,/x) = log A. (4.22) 



Moreover^ as N ^ oo^ the N -strip Gibbs measure (see Eqn con- 
verges weakly to a limiting probability distribution P that is represented by a 



17 



positive recurrent Markov chain with states (t,cr)^ the transition matrix 
P = {P((t,cr), (t',cr'))} and the invariant distribution tt = {7r(t,cr)} where 



P((t,cr),(t',a')) 



A'((t,o-),(t^o-o)(/)(t^o-o 

Acf>{t,a) 



with the norm 



2 

T-C 



7r(t,o-) = (/)(t,o-)(/)^(t,o-) /((/), </)^>^_^ 



4.2 



Proof of Theorem 



we have the convergence 



The spectral gap for K imphes that ^ ip E ^t-c^ 



hm -^K^^ = ((V',(/>)t-c)</> 

in the norm of space i^_c- Moreover, let 11 denote the operator of projection 
to the subspace spanned by the eigenvectors of K different from 0. Then 



1 



I nKP I 



T-C 



< 1 



hm 



(nKP) 



N 



= 0. 



T-C 



In turn, this imphes that 
1 



1 



^ log S7v(/x, /3) = ^ log trx-cK^ ^ log A. 



Convergence of the Gibbs measure follows as a corollary. □ 



5 Concluding remarks 

This paper makes a step towards determining the subcriticality domain for an 
Ising-type model coupled to two-dimensional causal dynamical triangulations 
(CDT). In doing so we employ transfer-matrix techniques and in particular 
the Krein-Rutman theorem. We complement the discussion of the previous 
sections with the following two concluding remarks: 

Remark 1. It is instructive to summarise the logical structure of the argu- 
ment establishing Lemma \4^ and Theorems \4.1\ and\4^' 
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Figure 3: Ag = A and A^ are the maximal eigenvalues of the matrix Q and 
a related matrix T respectively (see Appendix A). The area above the black 
curve is where the condition (4.20) holds true. 



• First, (4.21) holds iff condition (4.20) holds: see the proof of Lemma 



• Next, (4.21) implies that K Z5 a Hilhert-Schmidt operator and K is a 



trace class operator in 



The last fact, together with the property of positivity (4.16)^ allow us to 



use the Krein-Rutman theory, deriving all assertions of Theorems 14^1 
and \4^S\ 



On the other hand, if (4.20) fails (and therefore (4.21) fails), it does not 
necessarily mean that the assertions Theorems and 4-^ fail- I'^ other 
words, we do not claim that the boundary of the domain of parameters /3 



and ji where the model exhibits uncritical behavior is given by Eqn. (4.20). 
Moreover, Figure\B shows the result of a numerical calculation indicating that 



the condition (|4.20) is worse than (4.11) for (moderately) large values of /3. 



An apparent condition closer to necessity is the pair of inequalities (4.14) 
for some (possibly) large r. This issue needs a further study. 
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Remark 2. Physical considerations suggest that the critical curve in the 
{/3^g) quarter-plane would have some predictable patterns of behavior: as a 
function of it would decay and exhibit a first- order singularity at a unique 
point /3 = /3cr ^ (0, oc). 

A plausible conjecture is that the boundary of the critical domain coincides 
with the locus of points {(^^jj) where A looses either the property of positivity 
or the property of being a simple eigenvalue. This direction also requires 
further research. 
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A Proof of Lemma 14, IL 

By definition the trace ( |4.21[ ) we need to calculate the series 

tr(KTK) = J]KTK((t,o-),(t,o-)) 
(t,o-) 

= J] ir((t,^T),(t',0)^((t,^),(t',^')) 

(t,(T),(t',(T') 

= K\{t,a),{t\a')). (A.l) 

(t,<T),(t',<T') 

A single-strip triangulation t consists of up- and down-triangles. Accord- 
ingly, it is convenient to employ new labels for spins: if a triangle t{l) is an 
/th up-triangle then we denote it by t''^^] the corresponding spin a{j) will be 
denoted by cr(jp. Similarly, if t{j) is an /th down-triangle then we denote it 
by t^^; the spin a{j) will be denoted by ct^q. Consequently, the triangulation 
t and its supported spin-configuration cr are represented as 

t:=(tup,tdo) and <t := (<Tup, <Tdo)- 
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Here 
and 



(t , . . . , t )^ t^o (^do5 • • • 5 



^ up up5 • • • 5 up/ 

assuming that the supporting single-strip triangulation t contains n up- 
triangles and m down-triangles. (The actual order of up- and down-triangles 
and supported spins does not matter.) 



The same can be done for the pair (t^cr^) as illustrated in (A.l). Let 
recall that the triangulations t and are consistent (t ^ t') iff number of 
the down-triangles in t equals that of up-triangles in t^ 



m 



{t\a') 



{t,a) 



m ,1 


,2 m 




up ^ 


A'"^ 


A/ 






r / \ / 


^/\ 






up \ / 


^ up y 



,3 ,1 
^ up^ up 




Figure 4: Interaction betwen the elements z^, Tji 



^do^do 



To calculate the sum ( |A.l ) we divide the summation over (t^, a') into a 
summation over (t'^p,<p) and {t'^^^ cr'^o) - Firstly, fix a pair (t'^p,<p) and 
make the sum over (t^^^, cr^^^). Note that the term V{{t^ cr)^ (t^, cr')) depends 
only on CTdo and cr'^^. Consequently, 



^do'^ do 



(A.2) 



g-/3ff(0-)g-2/3y((t,<T),(t',<T'))g-/.n(t) ^ 



g-/3ff(0-')g-/.n(t')_ 
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The sum in the right-hand side of ( A.2[ ) can be represented in a matrix form. 
Denote by e±i the standard spin-1/2 unit vectors in M^: 



and 6-1 



Next, let us introduce a 2 x 2 matrix T where 
T = e-^ 



(A.3) 



Denote by i = 1, . . . , nup{t') the number of down-triangles in t' which 
are between the ith and {i + l)th up-triangles in t^ Let nup{t') = k then 



E 



g-/3if(O-0g-/in(t0 



^do'^ do 



where the matrix M is the sum of the geometric progression 



(A.4) 



n=l ^ ^ / 



(A.5) 



Using the same procedure we can obtain the sum over all up-triangles into 
the triangulation t. The only difference is the existence of marked up-triangle 
in the strip: let as before nup{t') = ndo{t) = k then 



k-i 



See Figure [4| for illustration of these calculations ( |A.4 ) and (A. 6). Further, 
supposing the existence of the matrix M and using (A.4) and (A. 6) we obtain 
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the following: 



^ g-/3H(<T)g-^,n(t) ^ ^-I3H((T') ^-nn(t') 



X 

*"P.<^up (t:io:0-do) 

^ g-2/3y((tdo,crdo),«p,cr'„p)) 



Necessary and sufficient condition for tlie convergence of the matrix se- 
ries for M is that the maximal eigenvalue of matrix T is less then 1. The 
eigenvalues of T are 

A± = e(^-'^)±e-(^+^), (A.8) 
and the above condition means that A+ < 1 or, equivalently, 

> ln(2cosh(/5)). (A.9) 



Under this condition (A.9), the matrix M is calculated explicitly: 



M = 



(1 _ e(/5-M))2 _ e-^n 



^e2/3 ^ (1 _ g4/3)g-(/3+p) 1 

^ I 1 e2/3 + (1 - e4/3)e-(''+'^) 



(A.IO) 



We are now in a position to calculate the sum in ( |A.l ). To this end, we 



again represent it through the product of transfer matrices. Pictorially, we 
express the above sum as the partition function of a one-dimensional Ising- 
type model where states are pairs of spins (a^^, aj^p) and the interaction is via 
the matrix T between the members of the pair and via matrix M between 
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neighboring pairs. More precisely, define the following 4x4 matrices: 



Q-- 



V 



Qr 



(2) 

++■'"++ 
(2) 

(2) 

+ '"'++ 

2/3^ ^(2) 



Qt 



Q 



tm 



( e2^t++m++ 
.(2) 



e~^^m m++ 

m m_+ 

(2) 

-2/3 (2) -2/3 (2) 

-2/3 (2) 



(2) 



-2/3 



m m 



(2) 



t_+m 

e-2/3i ^(2) 

e" 



m m 



(2) 



m m+_ 

e^^m m y 

(2) 

m+_m_i 

(2) 

e^^m m^ly 

t+_m 

t m+_ 

e^^t m j 



(A.ll) 



(2) 
++ 



(A.12) 



(A.13) 



(2) 



(2) 



-2/3^ _(2) 



t m^^\. e^l^t miy 



++ 



t+_ml 



(2) 



(A.14) 



where mij^rrr--^ and t^j {i^j G { — , +}) are elements of the matrices M, M^, 
and T respectively. 



(A.7) 



Now for the sum under consideration we obtain using representation 

g-2/3y((tdo,crdo),(t'up,cr;.p)) 



E 

(t,<T),(t',(T') 



A-^((t,<T),(t',<T')) 



X 



(tdo,crdo),(t;,p,cr;^p) 



.1=1 
k 



1=1 
k-1 



/=1 1=1 



1=1 
oo 



tr 



(A.15) 



k=0 



k=l 



By the construction the matrix Q is greater then Qt elementwise. Thus the 
eigenvalue of matrix Q is greater than the eigenvalue of the matrix Qt (it 
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follows from the Perron- Frobenius theorem). Therefore the necessary and 



sufficient condition for the convergence in (A.l) is that the largest eigenvalue 



of Q is less than 1. It is possible to calculate its eigenvalue analytically. In 



order to express the eigenvalues of Q it is convinient to use notations (4.18) 
and (4.19). In this notations the matrix M, i.e. (A. 10), is represented as 



following 



M 



m 
1 



1 



m 



The equations for the eigenvalues of Q are: 



Ai = 

A2 = 

As = 

A4 = 



c\^{m^ - 1) 
c^e-^im^ - 1) 



c2(m2 + 1) (cosh /3) 1- Wl- 



{m? — ly 



c^{w? + l)(cosh/3) {^^ ^ 



(cosh/3)2(m2 + l)2 

{m? - 1)2 
(cosh/5)2(m2 + 1)2 



A straightforward inspection confirms that the largest eigenvalue is given by 
A4. The condition A4 < 1 coincides with (A.l). This completes the proof of 
Lemma 14.11 



□ 
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